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8/4  Zhelyabov  str.,  Kyiv,  03057 
Ukraine 


Summary 


This  lecture  is  devoted  to  unsteady  processes  in  the  flows  with  natural  and  artificial  supercavitation. 
We  consider  two  classes  of  phenomena: 

1)  forced  non-stationarity  of  the  flow  which  is  induced  by  external  causes  - a model  velocity 
change,  an  ambient  pressure  impulse,  a variation  of  gas  supply  into  a cavity  etc.; 

2)  self-exited  oscillation  arising  due  to  internal  instability  of  gas-filled  supercavities. 

Results  of  computer  simulation  of  the  dynamic  supercavitation  processes  are  presented.  They  were 
obtained  by  using  the  approximate  mathematical  model  based  on  the  G.Logvinovich  independence 
principle  of  the  supercavity  expansion. 

A comparison  of  unsteady  behaviour  of  axisymmetric  and  two-dimensional  supercavities  is  given. 
Formulation  and  solution  of  two  2-D  problems  are  presented: 

1)  the  problem  on  instability  of  the  2-D  gas-filled  supercavity; 

2)  the  problem  on  evolution  of  the  2-D  supercavity  past  an  oscillating  wedge. 

The  solutions  are  based  on  the  M.Tulin's  linearized  cavitation  scheme. 

The  lecture  material  is  illustrated  by  sequences  of  motion-picture  frames  and  photographs  of  the 
unsteady  supercavitation  processes.  They  were  obtained  at  the  hydrodynamic  laboratory  of  the  Institute 
of  Hydromechanics  of  NAS  of  Ukraine  (IHM  UNAS). 


1 Main  types  of  dynamic  supercavitation  processes 


We  consider  main  types  of  the  supercavitation  flows  and  characteristic  for  each  of  them  unsteady 
phenomena  which  are  well  simulated  by  using  the  accepted  mathematical  model. 


1.1  SUPERCAVITY  FORMATION  DURING  HIGH-SPEED  WATER  ENTRY 


When  the  high-speed  motion  in  water  is  experimentally  investigated  at  the  IHM  UNAS,  models  are 
accelerated  with  the  vapor-gas  catapult  up  to  velocities  500  4-  1400  m/sec,  enter  into  water  through  the 
membrane  and  then  move  under  water  on  inertia  in  the  natural  super-cavitation  regime  [1,  2].  Length 
of  the  distance  is  35  m. 

Fig.  1 shows  consecutive  frames  of  the  initial  period  of  the  model  motion.  The  cell  of  the  coordinate 
mesh  is  0.2  x 0.2  m.  The  cavitator  diameter  of  model  is  Dn  =1.5  mm,  the  model  velocity  is  V0=  980 
m/sec.  One  can  see  that  the  cavity  formation  process  has  a stage  of  cavity  closure  caused  by  rapid 
increase  of  the  water  pressure  when  the  model  and  gases  penetrate  into  water.  In  this  case,  if  a stability 
of  the  model  motion  does  not  loss,  then  the  normal  formation  of  the  supercavity  goes  on  in  further. 


Paper  presented  at  the  RTOAVT  Lecture  Series  on  “Supercavitating  Flows  ”,  held  at  the  von  Karmctn 
Institute  (VKI)  in  Brussels,  Belgium,  12-16  February  2001,  and  published  in  RTO  EN-010. 
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1.2  Natural  high-speed  supercavities 


A main  characteristic  peculiarity  of  high-speed  super-cavities  is  their  very  huge  aspect  ratio  X = Lcl  Dc  = 
70  -s-  200,  where  Lc  is  the  cavity  length;  Dc  is  the  cavity  mid-section  diameter.  The  vapor  cavitation 
number  is  unique  similarity  parameter  for  such  flows: 


a = 2(^°°  ~fv) 

PV 2 


(1) 


where  is  the  pressure  at  infinity;  pv  = 2350  Pascal  is  the  saturated  water  vapor  pressure  (at  the 
temperature  20  °C);  p is  the  water  density.  This  type  of  flow  corresponds  to  the  cavitation  numbers  o < 
10  ~3 . 

Fig.  2 shows  a sequence  of  frames  of  shooting  the  high-speed  supercavity  {Dn  = 1.2  mm;  V = 1075 
m/sec;  Lc  = 18  m).  The  shooting  frequency  in  this  experiment  was  N=  4200  frames/sec,  the  time  interval 
between  the  frames  was  1 / N = 0.24  msec. 

A typical  unsteady  process  at  the  high-speed  motion  of  the  supercavitating  model  on  inertia  is  fast 
decrease  of  both  the  velocity  and  the  cavity  length.  The  super-cavitation  regime  of  motion  remains  until 
the  model  is  fully  placed  within  the  cavity. 

The  second  peculiarity  of  the  high-speed  super-cavities  is  the  cavity  closure  mechanism.  The 
cavities  in  the  described  experiments  fluently  closed  practically  in  a point  in  contrast  to  the  vapor 
cavities  when  cavitation  numbers  o <0.01.  For  latter,  the  unsteady  closure  with  periodic  arising 
reentrant  jet  is  characteristic.  It  is  seen  from  sequential  frames  of  the  experiments  that  the  wake  past  the 
cavities  has  periodic  bubble  structure.  Absolute  frequencies  of  the  periodic  structure  in  the  wake  were 
in  range  80  + 140  kHz  in  various  experiments  when  Dn  = 1.2  + 1.5.  In  Fig.  2,  the  frequency  is  equal  to 
140  kHz  . 

The  analysis  showed  that  one  of  causes  of  the  periodic  bubble  wake  formation  may  be  elastic 
vibrations  of  the  model. 

The  third  peculiarity  of  the  high-speed  supercavitation  motion  in  water  is  the  stability  mechanism  of 
motion.  For  the  supercavitating  models,  the  classical  condition  of  the  motion  stability  in  continuum 
is  not  fulfilled.  It  consists  in  location  of  the  point  of  the  hydrodynamic  force  application  past  the  center 
of  the  body  mass.  The  experiments  showed  that  the  stabilization  of  free  motion  of  the  supercavitating 
models  is  attained  owing  to  ricocheting  the  model  tail  from  the  inner  cavity  walls  when  V > 300  m/sec 
[1,2]. 

In  the  experiment,  action  of  this  stabilization  mechanism  appears  in  periodic  perturbations  of  the 
supercavity  surface  that  develop  according  to  the  independence  principle.  In  Fig.  3,  two  experimental 
photographs  are  showed:  the  model  at  the  time  of  touching  the  upper  internal  cavity  wall  (a)  and  the 
cavity  part  perturbed  owing  to  contact  with  the  model  (b).  Exposure  time  for  shooting  was  3 • 1 0 sec. 
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1.3.  Artificial  ventilated  cavities 


The  supercavitation  regime  may  be  created  in  water  at  the  moderate  velocities  V = 10  4-  100  m/sec.  The 
basic  similarity  parameters  of  the  ventilation  flow  are: 


pV2  ’ 


(2) 


where  pc  > pv  is  the  cavity  pressure;  Fr  is  the  Froude  number;  g is  the  gravity  acceleration;  Dn  is  the 

cavitator  diameter.  The  cavitation  numbers  1 0 “2  < a <0.1  correspond  to  this  type  of  flow. 

The  Froude  number  Fr  characterizes  the  distorting  effect  of  the  gravity  on  the  cavity  shape. 
Estimations  show  that  it  is  considerable  when  <r  Fr  < 2 [3], 

The  dynamic  similarity  parameter  yS  > 1 , that  is  equal  to  relation  of  the  vapor  cavitation  number  cry  and 
its  real  value  a , plays  an  important  role  at  calculation  of  the  unsteady  ventilation  flow. 

The  value  /3  = 1 corresponds  to  the  natural  supercavitation.  When  the  parameter  /3  increases, 
significance  of  elasticity  of  the  gas  filling  the  ventilated  cavity  increases  as  well. 


Fig.  3.  Cavity  perturbation  at  ricocheting: 
Dn  = 3 mm,  V = 690  m/  sec 


Fig.  4.  Ventilated  cavity  reaction  on  the  pressure  impulse: 
Dn  =20  mm,  V=  9 m/sec 


Fig.  4 shows  a sequence  of  motion-picture  frames  of  the  ventilated  cavity  deformation  under 
action  of  the  ambient  pressure  impulse.  The  pressure  impulse  in  water  was  created  with  the 
compressed-air  catapult.  The  air  bubble  created  by  the  catapult  is  seen  in  the  upper  part  of  the  frames. 
The  frames  show  that  the  pressure  impulse  results  in  the  axisymmetric  pinch  of  the  cavity. 

A characteristic  unsteady  process  for  ventilated  cavities  is  the  cavity  evolution  at  changing  or  ceasing 
the  gas-supply.  The  experiments  show  that  this  process  is  determined  by  type  of  the  gas-leakage  from  the 
cavity  [3,  4].  In  Fig.5,  the  typical  experimental  graphs  of  Lc(t),  Dc{t)  are  shown  when  the  air-supply 
was  instantly  turned  on  and  turned  off  [5].  The  portion  type  of  the  air-leakage  from  the  cavity  was 
maintained  during  this  test.  It  is  seen  that  the  cavity  length  increased  by  linear  law.  This  process  was 
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decelerated  only  at  approach  to  the  stationary  regime.  On  the  contrary,  the  cavity  decreased  with 
increasing  speed  when  the  cavity  closes. 
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Fig.  5.  Ventilated  cavity  history  at  air-supply  change 


1.4  Instability  of  ventilated  cavities 


Arising  the  self-induced  oscillations  of  the  ventilated  cavities  [6,  7]  is  interesting  unsteady  phenomenon. 
If  too  many  gas  is  supplied  into  the  cavity,  then  it  can  become  unstable.  In  this  case  waves  arise  on 
the  cavity  surface,  it  pulsates  along  the  its  length  and  width,  the  gas-leakage  from  the  cavity  is  realized 
by  separation  of  great  cavity  portions  (air  pockets).  A photograph  of  pulsating  axisymmetric  cavity  is 
presented  in  Fig.  6 (Yu.F.Zhuravlev). 

Theoretically,  this  phenomenon  was  explained  by  E.V.Paryshev  [8,  9]  on  the  basis  of  both  the 
G.V.Logvinovich  independence  principle  and  the  equation  of  the  mass  of  gas  in  the  cavity  balance. 


1.5  Cavities  after  the  water  entry  from  the  atmosphere 


The  air  cavity  forming  after  water  entry  of  bodies  from  the  atmosphere,  closes  on  the  depth  or  at  the 
water  surface  depending  on  the  initial  conditions  [10].  The  cavity  pressure  is  lower  than  atmospheric  one 
at  the  cavity  closure  time.  During  further  immersion  of  the  body,  the  cavity  fast  decreases  because  of 
both  the  static  pressure  growth  and  the  loss  of  the  air  entrapped  from  the  atmosphere. 

The  unsteady  process  of  the  depth  cavity  closure  is  well  described  by  the  approximated  mathematical 
model  based  on  the  independence  principle  [11]. 


Fig.  6.  Pulsating  axisymmetric  cavity 


We  have  investigated  experimentally  the  interesting  phenomenon  of  wave-shaped  deformation  of  the 
cavities  forming  after  vertical  water  entry  of  bodies  with  velocity  about  10  m/sec  [12]  (Fig.  7).  We 
showed  that  cause  of  this  phenomenon  consists  in  excitation  of  one  of  fundamental  frequencies  of  the 
cavity  filled  by  air  entrapped  from  the  atmosphere. 
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Fig.  7.  Cavities  after  water  entry  from  the  atmosphere 


2.  Mathematical  model  and  calculation  algorithm 


We  use  the  approximate  mathematical  model  based  on  the  G.V.Logvinovich  independence  principle  of 
the  cavity  section  expansion  for  computer  simulation  of  all  the  mentioned  unsteady  processes.  The 
mathematical  model  includes  the  following  equations: 

1)  The  equation  of  expansion  of  the  unsteady  cavity  section  [13]: 


d2S(r,t)  kxAp{x,t) 

dt2  p 


x (t)-l(t)<£<x(t), 


S(T,T) 


dS(z,r) 

dt 


(3) 


where  Ap{z,t)  = p„(£)  + P\(t)~  pc(t).  Here,  z <t  is  time  of  section  E,  creation;  x(t)  is  the  current 
absolute  .r-coordinate  of  the  cavitator;  /;,(7)  is  the  perturbation  of  the  ambient  pressure;  lit)  is  the 
cavity  length;  cx  is  the  cavitation  drag  coefficient;  kl  = An  I A2 ; A = 2 is  the  empirical  constant.  The 
hydrostatic  pressure  px  can  vary  from  one  section  to  another  when  the  body  moves  with  the  variable 
depth  in  ponderable  fluid. 

2)  The  equation  of  the  mass  of  gas  in  the  cavity  balance  when  the  gas  expands  by  the  isothermal  law 

[8]: 

^[pc(t)Q(t)]  = Plq,„-qout(t)l  (4) 

dt 


where  pc(t)  = pc(t)/<70 ; <r0  is  the  initial  cavitation  number;  Q is  the  cavity  volume;  qin  and  qou,(t)  are 
volumetric  rates  of  air-supply  to  the  cavity  and  air-leakage  from  the  cavity,  respectively,  referred  to  px . 
We  assume  that  the  cavity  pressure  pc{t)  changes  synchronously  along  the  cavity  length. 


3)  The  equation  of  rectilinear  motion  of  the  supercavitating  model: 


dU 

m 

dt 


(5) 
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where  m is  the  model  mass;  U(x,t)  is  the  model  speed;  XF  is  a sum  of  acting  forces  (cavitation  drag, 
propeller  thrust,  gravity  when  vertical  motion  etc.).  In  problems  on  the  supercavitating  model  dynamics, 
we  use  a complete  set  of  the  dynamic  equation  of  axisymmetric  model  [14]. 

This  mathematical  model  may  be  applied,  strictly  speaking,  when  the  parameters  change  not  very 
fast.  In  the  case  of  oscillation,  we  have  estimation  for  dimensionless  angular  frequency  k«2n  a/V„, 
where  a is  the  sound  speed  in  the  gas  filling  the  cavity.  It  follows  from  the  assumption  that  the  cavity 
pressure  pc{t)  changes  synchronously  along  the  cavity  length. 

We  have  developed  a "fast"  numerical  algorithm  and  a number  of  computer  codes  on  the  basis  of  the 
Eqs.  (3),  (4)  and  (5).  That  codes  allow  to  reproduce  the  unsteady  supercavitation  processes  of  the 
mentioned  types  on  a computer  screen  [2,  15,  16].  Examples  of  the  computer  simulation  are  given 
below. 
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Fig.  8.  Result  of  the  PCAV  code  operation:  cavity  reaction  on  the  water  pressure  impulse 


3.  Cavity  deformation  under  action  of  the  external  pressure  perturbation 


The  PCAV  code  (PerturbedCavity)  is  intended  for  computer  simulation  of  supercavity  deformation 
under  action  of  the  ambient  pressure  perturbations  in  water. 

Fig.  8 presents  the  computer  simulation  of  action  of  the  external  pressure  impulse  imitating  the 
underwater  explosion.  In  this  case  the  cavity  pinch  occurs  when  the  pressure  impulse  amplitude  is 
enough.  It  qualitatively  corresponds  to  the  experimental  shooting  in  Fig.  4. 

In  Fig.  9,  a result  of  computer  simulation  of  the  supercavity  formation  after  high-speed  water  entry 
through  the  solid  wall  is  shown.  The  extra  pressure  impulse  shape  pt(t)  in  Eq.  (3)  was  evaluated  by 
numerical  solving  the  problem  on  water  penetration  of  a solid  body  of  revolution.  In  each  time  its  shape 
coincides  with  the  cavity  shape  ("principle  of  freezing  the  free  boundaries")  . The  model  contour  is 
plotted  by  thin  lines,  the  cavity  contour  is  plotted  by  bold  lines.  The  cavity  contour  at  absence  of  the 
pressure  perturbation  is  shown  by  dashed  line  for  comparison.  In  this  case  the  pressure  impulse  arising 
at  the  water  entry  of  the  model  results  in  local  pinch  of  the  cavity  as  well.  A comparison  with  the  Fig.  1 
shows  qualitatively  agreement  of  the  calculation  and  experiment. 

The  same  kind  of  deformation  of  the  cavity  occurs  at  computer  simulation  of  water  entry  of  the  model 
from  within  an  expanding  gas  bubble  (Fig.  10).  In  this  case  the  pressure  in  the  fluid  may  be  computed 
exactly.  The  model  velocity  is  assumed  to  be  much  more  than  the  bubble  expansion  velocity.  The 
pressure  in  the  cavity  is  equal  to  the  pressure  in  the  bubble  for  the  starting  penetration  period. 
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: Current  x / Rn  = 291 .0  t V0 / Rn  = 291 .5 

x.'P.  240  160  80  0 

Fig.  9.  Result  of  PCAV  code  operation: 
high-speed  water  entry  through  a wall 


Fig.  10.  Result  of  PCA  V code  operation: 
water  entry  from  within  an  expanding  gas  bubble 


4.  Dynamic  properties  of  gas-filled  super-cavities 


We  describe  briefly  a procedure  of  the  stability  analysis  and  calculation  of  self-induced  oscillation  of  the 
axisymmetric  gas-filled  cavities  according  to  the  works  [8,  9,  14], 


We  make  dimensionless  the  Eq.  (3)  using  the  initial  cavity  length  l0  and  the  velocity  Vx  as  scales: 

d2S(r,t)  k^o(t) 


dt2 


t -l(t)<T  <t. 


(6) 


We  consider  for  simplicity  that  the  cavitator  is  negligibly  small  compared  to  the  cavity: 

S(t,t)  = S(t-l(t),t)->  0, 


(7) 


Twice  integrating  the  Eq.  (6)  with  taking  account  of  (7)  and  using  the  Dirichlet  formula  for 
interchanging  integrals,  we  obtain: 


k\  (70 


t-r-2 


i 

j"  (t  — u)o(u)du, 


(8) 


where  (J(t)  = a(t)/(J0 . Substituting  (8)  into  the  cavity  closure  condition  (7),  we  obtain  the  equation 
connecting  two  unknown  functions  (fit)  and  /(/): 

t 

l(t)  = 2 j \t — u)o(u)du  . (9) 

t-m 


The  dimensionless  equation  of  the  mass  of  gas  in  the  cavity  balance  (4)  has  the  form: 

dt 


(10) 


We  calculate  the  cavity  volume,  integrating  the  Eq.  (8)  along  the  cavity  length  and  using  again  the 
Dirichlet  formula: 


i\t) 


+ 


l 

J(t  — u)2a(u)du 


t-m 
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The  work  [3]  presents  the  empirical  law  of  the  air-supply  for  steady  axisymmetric  cavities  when  effect  of 
gravity  is  weak: 


</,«,= /KA 


'Z--\ 

v °o 


/ = 0.01  h-  0.02, 


where  Sc  = nlT  / 4 is  the  cavity  mid-section  area.  We  have  for  the  considered  processes  pv  « px  . For 
the  weakly  perturbed  unsteady  cavities,  we  accept  the  quasistationary  law  of  the  gas-leakage  of  the  same 
structure  [8]: 


4,nA0  = rscu) 


P 

o{t) 


(11) 


The  values  /)  = 1,  <7  = 1 correspond  to  the  natural  vapor  supercavity.  In  this  case  the  Eq.  (10)  is 
satisfied  identically. 


4.1  Instability  of  axisymmetric  ventilated  cavities 


A set  of  the  Eqs.  (9)  and  (10)  is  related  to  the  class  of  dynamic  system  with  distributed  lag  [17],  It  has  the 
only  stationary  point  <7  = l,  / = 1.  We  investigate  it  on  stability  relatively  to  small  oscillations. 
Representing  the  unknown  functions  in  the  form: 

dr(t)  = l + £al(t),  l(t)  = l+£l,(t),  e~o(  1), 

and  linearizing  the  equations  and  excepting  /,  (t) , we  obtain  the  uniform  equation  with  respect  to  <7,  (/) : 

i 

<7 1 ( t ) - 12(/3  - 1)  \d(9  - Dct,  (t  - 9)dd  + - jj8(2j8  - l)cr,  (t)  = 0.  (12) 

o 2 

Following  to  a usual  procedure  of  investigation  of  solutions  of  the  equations  on  stability  with  respect  to 
small  oscillations,  we  do  the  substitution  in  the  Eq.  ( 12) 

C71(t)  = a^',  p=X  + jk 

and  obtain  its  characteristic  equation: 


/t3+|^(2J8-l)Ju2+12(/3-l)[iu(^  + l)+2(e-/J-l)]  = 0.  (13) 

Here,  k is  the  reduced  oscillation  frequency,  k = (ol0  / ; X is  the  oscillation  increment.  The  oscillation 

with  frequency  k damps  when  X < 0 and  increases  indefinitely  when  X > 0. 
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Fig.  1 1.  Real  part  of  the 
characteristic  Eq.  ( 1 3)  roots 
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Fig.  12.  Modes  of  the  axisymmetric 
cavity  pulsation 


In  the  first  case  the  solution  of  the  Eq.  (12)  is  asymptotically  stable,  in  the  second  case  it  is  unstable. 
The  frequency  k which  corresponds  to  the  increment  value  X = 0 is  called  the  fundamental 
frequency  of  the  dynamic  system. 

The  Eq.  (13)  contains  two  physical  parameters  y and  p . When  y ^ 0,  it  has  the  double  root  u = 0 
which  transforms  into  threefold  one  when  y = 0.  It  is  not  difficult  to  make  sure  by  direct  checking  that 
the  corresponding  secular  solutions  [17]  do  not  satisfy  the  Eq.  (12).  When  y = 0,  the  Eq.  (13)  has  a 
series  of  the  pure  imaginary  roots,  i.e.  the  fundamental  frequencies: 


Fig.  1 1 presents  graphs  of  distribution  of  the  real  part  of  the  Eq.  (13)  roots  when  y = 0.  When  /3  — > , 

the  curves  asymptotically  approaches  to  the  a- axis.  We  suppose  that  the  frequency  corresponding  to  the 
maximal  linear  increment  X "survives"  among  severe  fundamental  frequencies  for  the  given  value  of  /3  . 
In  a result  we  obtain  estimation  of  the  mode  boundaries  of  the  cavity  self-induced  oscillation  (solid  lines 
in  Fig.  12). 


Thus,  when  y = 0,  the  cavity  is  asymptotically  stable  when  1 < /3  < 2.645  and  unstable  when  P > 
2.645.  The  value  y = 0 corresponds  to  the  case  when  the  mass  of  gas  in  the  cavity  is  constant.  An 
analysis  of  the  characteristic  Eq.  (13)  shows  that  it  is  possible  to  point  a finite  interval  of  changing  the 
value  /3  for  each  0 < y < 0.08.  Outside  this  interval,  the  zero  solution  of  the  Eq.  (12)  is  asymptotically 
stable,  and  within  this  interval  it  is  unstable  (Fig.  13).  Within  the  instability  zone  and  when  y^  0,  the 
characteristic  equation  has  the  finite  number  of  pure  imaginary  roots.  Their  quantity  decreases  with 
increasing  y . When  y > 0.08,  roots  are  absent  in  general,  i.e.  the  cavity  is  asymptotically  stable  for  any 
P. 


The  established  behaviour  of  effect  of  the  gas-leakage  variability  on  the  cavity  stability  qualitatively 
explains  the  experimental  fact  that  the  cavity  pulsation  ceases  when  the  gas-supply  becomes  very  great. 
So,  in  the  experiments  [6,  7],  five  waves  on  the  pulsating  cavities  (N=  5)  have  been  observed  at  most. 
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Fig.  14.  Result  of  the  PULSE  code  operation  (mode  II) 


4.2  Self-induced  oscillation  of  ventilated  cavities 

For  the  first  time,  calculation  of  self-induced  oscillations  developing  after  the  stability  loss  of  the 
ventilated  cavity  was  performed  by  E.V.Paryshev  [9],  We  have  developed  the  computer  code  PULSE  that 
allows  to  reproduce  the  cavity  pulsation  on  a computer  screen  "in  real  run-time"  and  to  carry  out  its 
Fourier  analysis  [15]. 

Calculation  showed  that  pulsation  arises  in  the  dynamic  system  described  by  Eqs.  (3)  and  (4)  when 
the  parameter /3  is  in  the  linear  instability  zone  ( jS  > /3,=  2.645)  . The  oscillation  develops  until  the 
establishment  of  the  periodic  or  the  quasiperiodic  pulsation  with  discontinuous  dependence  l(t).  Fig.  14 
shows  a view  of  the  computer  screen  during  run-time  of  the  code  PULSE.  Fig.  15  shows  evolution  of 
spectrum  of  the  cavity  pressure  when  the  bifurcation  parameter  q0  = f50q in  increases.  Increasing  the 
bifurcation  parameter  q0  is  accompanied  by  spasmodic  appearance  of  new  frequencies  and  their  linear 
combinations  in  the  spectrum  pc . This  corresponds  to  passage  up  to  higher  modes  of  the  cavity 
pulsation.  The  spectrum  pc  becomes  complicated  within  each  mode,  remaining  the  line  one.  In  this  case 
the  basic  harmonic  changes  weakly. 

The  same  behaviour  of  unsteady  ventilated  cavities  is  observed  experimentally  [6,  7]. 

Fig.  16  shows  spectra  of  the  cavity  pressure  oscillation  under  action  of  forced  oscillation  of  the 
external  pressure  pi(t)  = Ksinkft  ( q(1  = 0.8,  K=  0.1).  The  corresponding  spectrum  pc(t)  when  the 

external  perturbation  is  absent  is  given  in  Fig.  15,  a.  We  have  the  modulation  when  kf  «k,  (a),  the 

synchronization  when  kf  ~k,  (b)  or  the  "chaotization"  when  k,  > k,  (c)  of  the  periodic  mode  1 in 

depending  on  a ratio  of  the  forcing  frequency  kt  to  the  cavity  oscillation  frequency  kl  (they  both  are 

referred  to  the  average  cavity  length  lm ). 
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Fig.  15.  Power  spectral  density  of  pL,(t)  at  the  self-induced  oscillation 
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Fig.  16.  Power  spectral  density  of  pcif)  at  the  forced  oscillation 


4.3  Wave  formation  on  the  cavity  after  water  entry  from  the  atmosphere 

We  already  spoke  of  the  wave  formation  phenomenon  on  the  cavities  forming  at  vertical  water  entry 
of  bodies  through  the  free  surface  from  the  atmosphere  with  the  velocities  5-^10  m/sec  [12]. 

It  was  established  experimentally  [18]  that  the  free  water  surface  influence  on  the  drag  coefficient  cx 

and  on  the  cavity  shape  near  the  cavitator  propagates  only  on  the  depth  of  1.5  + 2 cavitator  diameters. 
This  result  gives  background  to  use  the  approximate  Eq.  (3)  for  computer  simulation  of  penetration  of  the 
bodies  through  the  free  surface  into  water. 

To  explain  this  phenomenon  we  have  applied  results  of  the  linear  theory  of  gas-filled  supercavity 
instability.  We  showed  that  its  cause  consists  in  excitation  of  one  of  the  fundamental  frequencies  of  the 
cavity  filled  by  air  entrapped  from  the  atmosphere.  In  the  case  of  enough  high  both  the  body  mass  and  the 
initial  Froude  number  Fr0 , we  obtain  a priori  estimation  of  the  wave  number  N on  the  cavity: 


, B = (15) 

where  Eu0  — 2palm  / pF02  is  the  Euler  number.  The  estimation  (15)  is  in  good  agreement  with  the 
experimental  data. 

The  calculation  model  of  that  process  after  the  cavity  depth  closure  includes  three  Eqs.  (3),  (4)  and 
(5).  Their  solution  is  sought  numerically  by  iteration  process  for  each  time  step.  Fig.  17  presents  graphs 
of  both  the  cavity  length  and  the  cavity  pressure  as  functions  of  the  cavitator  immersion  x.  The  graph  of 
the  dependence  / (x)  which  was  calculated  without  taking  account  of  the  air  elasticity  is  plotted  by  a 
dashed  line.  A comparison  of  the  calculated  shape  and  the  experimental  shape  of  the  cavity  for  three 
sequential  times  is  given  in  Fig.  18. 
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The  obtained  good  agreement  of  the  calculation  and  experiment  for  such  complicated  unsteady  process 
confirms  adequacy  of  the  accepted  approximate  mathematical  model  of  the  unsteady  axially  symmetric 
supercavity. 


Fig.  17.  Cavity  history  after  vertical  water  entry:  graphs  of  pc  (t)  and  1(f) 


5.  Control  of  the  ventilated  cavities 


A problem  of  control  of  the  ventilated  cavity  is  the  problem  of  maintaining  the  cavity  dimensions  or  of 
varying  the  cavity  dimensions  according  to  the  given  law  by  regulating  the  gas-supply  into  the  cavity. 
Examples  of  the  typical  control  problems  are: 

1)  How  should  the  gas-supply  rate  be  varied  for  maintenance  of  the  invariable  cavity  if  the  motion 
velocity  and/or  the  motion  depth  change? 

2)  What  is  the  law  governing  the  cavity  collapse  if  the  gas-supply  to  the  cavity  is  rapidly  stopped? 

The  difficulty  of  the  problem  on  the  ventilated  cavity  control  is  caused  by  non-linearity  and  non- 
monotonicity of  the  dependence  qin(t)  = f (a)  [3,  4]  and  by  the  multiparametric  nature  and  lack  of 
knowledge  of  the  gas-leakage  laws  as  well. 


Fig.  18.  Cavity  history  after  vertical  water  entry: calculation, experiment 
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We  have  developed  the  computer  code  ACAV  (ArtificialCavitation)  for  computer  simulation  of  unsteady 
processes  at  the  ventilated  cavity  control.  It  allows  to  use  any  laws  of  varying  the  air-supply,  the  air- 
leakage  and  the  model  velocity  as  well. 

Fig.  19  shows  a view  of  the  computer  screen  during  run-time  of  the  code  ACAV.  In  this  case,  the 
model  quasi-stationary  law  of  the  gas-leakage  from  the  cavity  was  used: 

• 2 • 10-5 

?«*(0  = — J—  • (16) 

o\t) 

It  corresponds  to  the  second  type  of  the  gas-leakage  by  vortex  tubes. 

in  the  code  ACAV,  the  set  of  the  Eqs.  (3)  and  (4)  is  integrated  with  the  constant  step  with  respect  to 
dimensionless  longitudinal  variable  x.  During  computation  in  Fig.  19,  the  gas-supply  rate  into  the  cavity 
increases  in  2.4  times  on  the  distance  Ax=  0.5,  then  remains  constant  on  the  distance  Ax  =1.0  and  than 
decreases  to  its  initial  value.  The  calculation  gives  a characteristic  lagging  reaction  of  the  cavity  on  the 
changing  the  rate  which  is  observed  in  experiments. 


Fig.  19.  Result  of  the  ACAV  code  operation 


5.1  Cavity  reaction  on  changing  the  gas-supply  rate 


Firstly  we  consider  peculiarities  of  the  ventilated  cavity  reaction  on  changing  the  air-supply  rate 
qin  (t)  when  the  model  velocity  is  constant:  V0  = const.  Here  we  present  results  of  computer  simulation 
for  two  types  of  changing  the  supply:  1)  impulse  increase  of  the  air-supply;  2)  rapid  single  decrease  of 
the  air-supply. 

At  the  calculation,  we  used  the  following  parameters  of  the  model:  the  disk  cavitator  with  diameter 
Dn  = 200  mm;  the  cavitator  slope  angle  <5  = 10°;  the  model  tail  diameter  Ds=  350  mm;  the  model 
length  L = 7 m.  Also,  we  used  the  following  common  starting  values  of  the  parameters:  V0  =100  m/c; 
the  motion  depth  H = 5 m;  <70  = 0.04;  the  model  pitch  angle  \ff=  2°.  In  this  case  the  starting  cavity 
length  is  Lc  = 9.06  m when  the  balanced  air-supply  coefficient  is  q0  =7.813. 

In  Figs.  20,  21,  the  following  legends  are  used  for  marking  the  curves:  1 - q(t) ; 2 - a(t) ; 3 - /(/); 
4-  V(t). 

Fig.  20,a  demonstrates  a calculation  result  of  the  cavity  reaction  on  the  impulse  increase  of  the  air- 
supply.  One  can  see  that  the  cavity  length  reacts  on  the  rate  change  with  greater  lagging  than  the  cavity 
pressure.  All  the  cavity  parameters  return  to  their  starting  values  after  the  impulse  end. 
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A calculation  result  of  the  cavity  reaction  on  the  rapid  decreasing  the  air-supply  is  shown  in  Fig.  20,  b. 
In  this  case  the  cavity  does  a number  damped  oscillations  (a  transient  process),  then  establishes  new 
balanced  values  of  the  parameters:  Lc  =6.09  m;  o =0.06  m;  qin  = 1.563. 

A cause  of  the  transient  process  is  the  elasticity  effect  of  the  gas  filling  the  cavity  when  the  cavity 
volume  rapidly  decreases.  We  note  that  self-induced  oscillations  of  the  cavity  may  arise  for  other  starting 
parameters  (smaller  values  of  V0  and  higher  values  of<70 ). 


5.2  Cavity  reaction  on  changing  the  model  velocity 


Now  we  consider  peculiarities  of  the  ventilated  cavity  reaction  on  varying  the  model  velocity  V when  the 
gas-supply  rate  is  constant:  qin  = const. 

The  calculation  result  of  the  ventilated  cavity  reaction  on  increasing  the  model  velocity  is  shown  in 
Fig.  21,  a.  One  can  see  that  the  cavity  length  increases  very  significantly  (in  quasi-stationary 
approximation  it  is  proportional  to  the  velocity  square),  but  the  cavity  pressure  changes  weakly.  It  is 
typical  for  this  case  that  the  unsteady  cavity  length  considerably  exceeds  the  new  balanced  level  Lc  = 
13.0  m. 

The  calculation  result  of  the  ventilated  cavity  reaction  on  decreasing  the  model  velocity  is  shown  in 
Fig.  20,  b.  A comparison  with  Fig.  20,  a shows  a difference  from  the  cavity  reaction  on  the  air-supply 
decrease  that  consists  in  absence  of  the  transient  process.  The  new  balanced  cavity  length  is  Ic=2.34m. 

We  can  investigate  with  the  ACAV  code  the  ventilated  supercavity  behavior  when  the  model  motion 
velocity  V[t)  and  the  air-rate  into  the  cavity  qin  (/ ) change  simultaneous.  In  this  case  we  can  try  to 
obtain  answers  for  two  practically  important  questions: 

1)  May  we  accelerate  the  cavity  development  on  the  accelerating  phase  of  the  body  motion  by  means 
of  increasing  the  gas-supply  into  the  cavity  ? 

2)  Flow  long  may  we  maintain  the  cavity  with  given  dimension  by  means  of  increasing  the  gas- 
supply  when  the  body  velocity  decreases? 

The  computer  simulation  allows  to  conclude  that  the  air-supply  increase  when  the  model  velocity 
increase  does  not  result  in  considerable  decrease  of  time  of  the  cavity  development  to  the  balanced  length 
at  the  accepted  law  of  the  gas-leakage  from  the  cavity  ( 16).  This  conclusion  is  valid  in  the  case  when  the 
model  velocity  is  already  high,  and  the  supercavity  already  exists. 


Fig.  20.  Cavity  reaction  on  varying  the  gas-supply  rate 
a - impulse  increasing  the  gas-supply  rate;  b - rapid  decreasing  the  gas-supply  rate 
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Fig.  21.  Cavity  reaction  on  varying  the  model  velocity:  a - increasing  the  velocity; 

b - decreasing  the  velocity 


5.3  Influence  of  the  parameter  p and  the  gas-leakage  type 

Considerable  difference  between  the  natural  vapor  and  ventilated  supercavities  appears  in  the  case  of 
the  unsteady  flow.  In  this  case  the  parameter  p is  of  importance  that  characterizes  effect  of  elasticity  of 
the  gas  filling  the  cavity. 

A dependence  of  changing  the  cavity  length  Lc  on  the  parameter  /30  is  shown  in  Fig.  22.  The 
calculation  was  performed  with  the  code  ACAV  when  Dn  = 200  mm.  The  free  cavity  closure  and  the 
gas-leakage  law  (16)  were  accepted  at  calculation.  The  cavitation  number  was  constant  and  equal  to  a = 
0.06.  The  starting  cavity  dimensions  were:  Lc  = 6.04  m,  Dc  = 0.74  m. 

The  model  velocity  firstly  increases  on  50%  along  the  distance  10  m,  then  remains  constant  along  the 
distance  100  m,  next  decreases  to  its  starting  value.  A graph  of  the  model  velocity  variation  is  shown  by 
dotted  line. 

The  curve  1 corresponds  to  the  natural  cavitation  regime  /3  = 1.  In  this  case  the  starting  model 
velocity  was  VQ  =69.5  m/sec. 

In  the  artificial  cavitation  regime,  changing  the  parameter  ft  at  the  constant  a is  attained  by 
changing  the  starting  model  velocity  V0 . We  used  the  data: 

curve  2 - V0  = 50  m/sec,  P0  = 1.96; 

curve  3 - V0  = 30  m/sec,  PQ  =5.45. 

The  presented  graphs  visually  demonstrate  increasing  the  effect  of  the  gas  elasticity  with  increasing  the 
parameter  [i  . In  the  last  case  (the  curve  3)  undamped  self-induced  oscillations  of  the  cavity  develop 
according  to  the  theory  of  instability  of  gas-filled  cavities.  In  this  case  the  starting  value  of  the  parameter 
is  P = 7.358  > 2.645,  i.e.  the  cavity  is  unstable  according  to  (14). 
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Fig.  22.  Influence  of  the  parameter  ji 
the  cavity  behaviour 


Fig.  23.  Cavity  behaviour  at  various  on 
types  of  the  gas-leakage 


In  Fig.  23,  the  dependence  Lc(t)  which  has  been  calculated  for  different  laws  of  the  gas-leakage  from 
the  cavity  is  shown.  The  model  velocity  changed  in  the  same  way  as  in  Fig.  22.  In  calculation  we 
assumed  that  the  leakage  coefficient  is  a single-valued  function  of  the  cavity  pressure:  Q = Q (pc ) . 
Different  curves  in  the  figure  correspond  to  the  following  gas-leakage  laws: 

curve  1 - the  1st  type  of  the  portion  gas-leakage  ; 

curve  2 - the  2nd  type  of  the  gas-leakage  by  vortex  tubes,  the  formula  (16). 

For  the  1st  type  of  the  gas-leakage  we  used  a function  obtained  by  the  experimental  data  [4] 
approximation: 


0 = 


5(7 

100(72  + f 


(17) 


Considerable  difference  of  the  cavity  behaviour  confinns  importance  of  observance  of  the  gas-leakage 
type  similarity  at  both  the  physical  modeling  and  the  computer  simulation  of  unsteady  ventilated  cavities. 

Presence  of  the  model  body  within  the  cavity  decreases  active  volume  of  the  cavity  filled  by  gas.  This 
does  not  influence  on  the  cavity  dimension  in  the  case  of  the  steady  flow  at  the  cavity  free  closure  (if  we 
do  not  take  into  account  the  possible  interaction  between  the  model  and  the  reentrant  jet).  In  the  non- 
stationary case  we  have  quite  another  situation.  It  is  necessary  to  take  into  account  decreasing  the  active 
cavity  volume  in  the  equation  of  the  mass  of  gas  balance  (4).  This  causes  to  changing  the  dynamic 
cavity  behaviour. 


6.  Comparison  of  unsteady  behaviour  of  2-D  and  axisymmetric  supercavities 


It  follows  from  said  before  that  the  approximate  equation  of  the  cavity  section  expansion  (3)  based  on  the 
G.V.Logvinovich  independence  principle  gives  an  effective  (and  usually  unique)  method  of  calculation 
of  the  unsteady  supercavitation  flows. 

The  independence  principle  is  approximate  one,  and  its  verification  has  great  methodological 
significance.  The  experimental  tests  justify  its  validity  in  wide  range  of  the  parameters  [2]. 

The  Eq.  (3)  well  describes  the  mean  part  of  the  supercavity,  i.e.  it  is  in  essence  linear  one.  A 
comparison  of  calculation  results  of  the  unsteady  cavity  shape  by  the  Eq.  (3)  with  the  theoretical 
solutions,  which  are  obtained  in  the  linearized  theory  of  two-dimensional  super-cavitation  flows  [19-22], 
is  of  interest. 
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We  use  the  solutions  of  two  2-D  problems  for  comparison: 

1)  the  problem  on  instability  of  the  2-D  gas-filled  supercavity; 

2)  the  problem  on  evolution  of  the  2-D  supercavity  at  the  forced  oscillations. 

The  linear  theory  application  to  this  problems  is  correct,  because  the  ventilated  cavities  always  have 
very  great  aspect  ratio  at  the  time  of  the  stability  loss  [6,  7], 

6.1  Mathematical  models  of  the  2-D  unsteady  supercavity 


Two  alternative  approaches  are  applied  to  solve  the  steady  and  non-stationary  2-D  linearized  problems  on 
the  supercavitation  flow  [23,  24]: 

1)  the  method  of  boundary  value  problems  for  analytical  functions; 

2)  the  method  of  integral  equations  (a  variant  of  the  method  of  boundary  integral  equations). 

The  first  method  gives  a solution  in  the  form  of  quadratures.  It  is  convenient  for  analytical 
investigation  of  stability  of  the  2-D  gas-filled  cavity  . 

The  second  method  is  more  convenient  for  numerical  calculations  when  the  cavity  length  is  variable. 
It  is  easy  generalized  for  the  bounded  flows,  for  the  super-cavitating  hydrofoils  in  hydrodynamic 
cascades  and  for  the  hydrofoils  of  finite  span  as  well. 


6.1.1  Method  of  boundary  value  problem  for  analytical  functions 


We  construct  firstly  a solution  of  the  linear  unsteady  problem  by  the  1st  M. Tulin's  scheme  (Fig.  24).  For 
simplicity,  we  consider  so  called  "pure"  supercavity,  i.e.  we  assume  that  the  cavitator  dimension  is 
negligibly  small  compared  to  the  cavity. 
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Fig.  24.  Two-dimensional  "pure"  unsteady  supercavity 
by  the  1st  M. Tulin's  scheme 


All  the  variables  are  assumed  to  be  dimensionless.  The  initial  cavity  length  /0  and  the  velocity  are 
accepted  as  the  length  scale  and  the  velocity  scale,  respectively.  We  consider  two  analytical  functions: 
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the  complex  velocity  V(z,t)=(px-i(p  and  the  complex  acceleration  potential  &(z,t)=d+i%.  We 
have: 


p at  ax 

On  the  x-axis,  the  components  of  &(zj)  and  V ( z,t ) are  connected  by  relations: 

X 

Q{x,t)=9{-c° , t)+  (px(x,t)+  j(pxt(s,t)ds. 


(18) 


(19) 


X 

X(x,t)=  x(-°° , t)-(py  (x,t)~  J (pyt  ( s,t)ds . 


(20) 


The  cavity  pressure  pc(t)  must  be  constant  at  each  time  along  the  cavity  length.  Thus,  we  have  the 
boundary  condition  for  the  acceleration  potential: 

e(x,f)  = — , 0<j c </(*).  (21) 

The  cavity  length  l(t)  is  an  unknown  time  function.  The  cavity  pressure  pc  and,  therefore,  the  cavitation 
number  a are  unknown  time  functions  in  the  case  of  the  gas-filled  supercavity. 

Differentiating  the  Eq.  (2 1 ) with  respect  to  x and  applying  the  inverse  operator  N~] , we  obtain  the 
boundary  condition  for  perturbed  horizontal  velocity: 

(px(x,t)=C(t-x),  0 <x<l(t), 

where  C'(f)  is  an  arbitrary  time  function.  Out  the  interval  [0,  /(f)]  we  have  tpr  =0  because  of  the  flow 
symmetry  and  absence  of  sources  in  the  cavity  wake. 

According  to  the  1st  M.Tulin  scheme,  the  "pure"  super-cavity  has  a shape  of  ellipse  in  the  stationary 
case: 


F(x)U t0  =-^y/x(/0-x),  0 <x<l0. 

In  the  points  z = 0,  z = l0 , the  complex  velocity  has  singularities  of  the  order  of  x 

We  assume  that  the  solution  must  have  singularities  of  the  same  order  for  unsteady  supercavity  in  each 
time  on  the  ends  of  the  interval  [0,  /(f)] . 

Thus,  we  have  the  mixed  boundary  value  problem  for  the  function  V (zj)  in  the  upper  half-plane 
Imz  >0 : 


(pv=0,  at  -oo<x<0,  (px=C(t\  at  0<x<l(t),  (pr=  0,  at  l{t)<x<°°.  (22) 

The  solution  of  this  problem  in  the  class  of  functions  unlimited  on  the  ends  of  interval  [0,  /(f)]  is  given  by 
Keldysh-Sedov  formula  [25]: 


K(z,f) 


z(z  - 1) 


+ 


7T  V Z-l  n 


l -s  ds 


s — z 


(23) 


Here  and  later,  functions  V and  C are  related  to  cr0 . 

The  solution  (23)  gives  a logarithmic  singularity  of  the  pressure  at  infinity.  It  is  caused  by  variability 
of  the  2-D  closed  cavity  in  unbounded  fluid.  This  methodological  paradox  disappears  when  the  free 
boundary  is  in  the  flow. 
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To  determine  the  cavity  shape  F(x,t)  we  use  the  linearized  kinematic  condition  on  the  cavity 
boundary  [24] 

(pv(x,t)=  N F(x,t\  0 <x<l(t).  (24) 

From  here,  using  the  relation  (20)  and  the  Dirichlet  formula  for  interchanging  integrals,  we  obtain  the 
expression  for  both  the  shape  and  the  volume  of  the  unsteady  symmetric  cavity: 


.1  X 

F{x,  t)  = | (x  - s)xt  ( s,t-x  + s)ds  - ^x(sF  ~ x + s)ds, 
o o 

/ 

Q(t)=  — 2 J (/  - s)x(s,t  -l  + s)ds. 


(25) 

(26) 


where 


X{x,t)=~ 


[ a q 2 
■ + — + — In 


q(t)- 2n 


4\x(/-x)  ^2  n J] 

77  1 


x q 

arctan 

/ - x n 


1 - x 


l-s 


ds 


It  is  necessary  to  fulfill  the  natural  boundary  condition  (21)  to  eliminate  the  auxiliary  function  C(t). 
Calculating  the  acceleration  potential  on  the  interval  [0,  /(?)]  by  the  formula  (19),  we  obtain  the  equation: 


/ 1 - / r C(t  - s)ds 


■J 


■iln4=  = 2. 


(27) 


477  n J 7,v(/  - s)  n 2 

The  constructed  solution  of  the  unsteady  problem  (23),  (27)  contains  two  unknown  time  functions  ail) 
and  lit).  To  determinate  them  we  use  the  condition  of  the  cavity  to  be  close  in  each  time: 

/ / 

J (/  - s)x,  (s,t  - 1 + s)ds  - J x(sF  ~ l + s)ds  = 0 (28) 

0 0 

and  equation  of  the  mass  of  gas  in  the  cavity  balance  (4). 

In  the  case  of  steady  flow,  the  solution  (23)  has  the  form: 


V{z)= x 


/« 


2yjz(z-l0) 


+ 1- 


z-L 


(29) 


We  note  that  the  parameters  aQ  and  lQ  are  independent  for  stationary  "pure"  supercavity. 


Fig.  25.  Unsteady  flow  around  supercavitating  wedge 
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6.1.2.  Method  of  integral  equations 


Now  we  construct  the  solution  by  the  method  of  integral  equations  for  the  unsteady  problem  on 
supercavitation  flow  around  a thin  symmetric  wedge  with  unit  length  (Fig.  25).  Let  the  unsteady 
perturbation  of  the  flow  caused  by  deformation  of  the  wedge  sides  so  that  a shape  of  the  wedge  is  given 
by  the  function: 

y = ±f(x,t)~0(e),  0 < x < 1,  (30) 

where  e is  the  small  parameter.  In  this  case  the  kinematic  boundary  condition  on  the  upper  side  of  the 
wedge  is  added  to  the  boundary  conditions  (22),  (24): 

(pv  = N f(x,t),  0 < x < 1;  y = + 0.  (31) 

Because  of  the  flow  symmetry,  we  can  obtain  the  flow  satisfying  the  boundary  conditions  (24),  (31) 
when  we  place  a layer  of  2-D  sources  with  the  linear  intensity  q(x,t)  on  the  interval  [0,  /(f)]  of  the  x-axis. 
At  the  arbitrary  point  (x,  y)  of  the  flow,  the  sources  induce  the  total  potential: 


‘V) 

(p(x,y,t)=—  jq(s,t)lnj(x-sf  +y2  ds. 


The  source  intensity  is  equal  to  a jump  of  the  normal  velocity  of  the  fluid  cpv  at  passage  through  the  x- 
axis.  We  have  for  the  flow  symmetric  about  the  x-axis: 

q(s,t)  = 2(py(s,t),  0 <s<l(t).  (33) 

Applying  the  linear  operator  N to  the  expression  (32),  we  obtain  the  acceleration  potential  of  unsteady 
sources  distributed  along  an  interval  with  a variable  length: 


j m 

0{x,y,t)  = — J q(s,t) 


( x-sY+y 2 In  dt 


1 d I 

ds- 1 J q(s,t)\n^(x  - s)2  + y2  ds. 


Passing  to  the  limit  y — » 0 in  Eq.  (34)  and  substituting  it  in  the  boundary  condition  (24),  we  obtain  the 
integro-differential  equation 


40  , , , 40 

f q(sj)ds  | 3 f 

J x — .S'  dt  J 
0 0 


(t) 

J q(s,  t ) In  | x - 5 1 ds  -na(t)  = 0 , 


where  1 <x  < l(t).  The  intensity  of  the  sources  distributed  along  the  interval  [0,  1]  is  known: 


q(x,t)  = 2Nf(x,t),  0<x<l. 

Thus,  the  Eq.  (35)  may  be  finally  rewritten  in  the  fonn: 

/(P  q(s,t)ds  5 /<r  ) | , 

j H—  J q(s,t)\n  | x - s | ds  -no(t)  = Al  (x,f),  (36) 

j x — 5 dt  ] 

where  the  function  Ax  (x,  t)  in  the  right  part  is  easy  calculated  for  the  concrete  law  of  unsteady 
deformations  of  the  wedge  (30): 


• Nq(s,  t)ds 


t 

N q(s,t) hi  | x-s  | ds. 
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Then,  we  obtain  the  equation  of  the  upper  cavity  boundary  from  the  kinematic  condition  on 
the  cavity  (24)  when  1 < x < l(t): 

1 ^ 

F(x,t)  = N~l(p  v(x,t)  = — [q(s,t  — x + s)ds.  (38) 

2o 

Since  the  cavity  length  l(t ) and  the  cavitation  number  a(i)  are  unknown  time  functions  in 
the  general  case,  it  is  necessary  to  add  two  relations  to  the  Eq.  (36): 

1)  the  condition  of  solvability  of  the  Newmann's  external  boundary  value  problem  for  the 
velocity  potential 

m 

f q(s,t)ds  = 0;  (39) 

o 

2)  the  equation  of  the  mass  of  gas  in  the  cavity  balance  (4)  for  gas-filled  cavities.  In  the  case 
of  natural  vapor  supercavities  the  Eq.  (4)  is  replaced  by  the  condition  pc  = const. 

As  a result,  we  obtain  a set  of  three  equations  to  determine  the  functions  q(x,t),  l(t)  and  a(t) . 
It  must  be  integrated  with  respect  to  time  with  the  initial  conditions 

q(x,0)  = q0(x),  1(0)  = l0,  a(  0)  = <r0.  (40) 


The  condition  (39)  ensures  limitation  of  the  pressure  at  infinity.  In  the  partial  case  of  the 
steady  flow,  the  Eq.  (39)  is  the  cavity  closure  condition.  In  the  case  of  the  unsteady  flow,  the 
cavity  is  unclosed. 

It  is  necessary  to  distribute  2-D  vortexes  with  the  linear  intensity  y(x, t)  along  the  projection  of 
the  foil  on  the  x-axis  [24]  together  with  the  sources  in  the  general  case  of  nonsymmetric  flow 
around  the  supercavitating  hydrofoil.  This  adds  one  more  singular  integral  equation  to  the  set  of 
equations.  Our  paper  [26]  presented  results  of  solving  the  problem  on  supercavitating  flow 
around  an  oscillating  hydrofoil  near  the  water  surface. 


6.2.  Instability  of  the  2-D  ventilated  cavity 

Using  the  solution  of  the  equation  set  (27),  (28)  and  (4),  we  investigate  stability  of  the  "pure"  2-D 
supercavity  in  the  unbounded  flow. 

According  to  the  experimental  data  [6]  for  the  steady  plane  supercavities,  the  air-supply  rate  may  be 
approximated  by  the  linear  function  (in  the  dimensional  form): 

4,„  =Yb()vi\-  — \ (41) 

where  y is  the  empirical  coefficient;  b0  is  the  cavity  mid-section.  It  is  valid  for  the  steady  "pure"  cavity 
(29)  that 

b0= 

0 2 
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We  assume  that  when  oscillation  of  the  cavity  pressure  is  small,  the  air-leakage  rate  depends  in  quasi- 
stationary way  on  pc  (?) : 


Qout 


(l)=Yb(t)V„ 


1- 


g(0 

(T„ 


(42) 


Passing  to  the  dimensionless  variables  and  using  the  quasi-stationary  dependence  b(c r),  we  obtain  from 


the  Eq.  (4): 


at  2 


1-1-1 

P ° 


— YI 
1-1 
P 


(43) 


The  set  of  equations  (27),  (28)  and  (43)  relates  to  the  class  of  nonlinear  autonomous  equations  with 
distributed  lag  of  Volterra  type.  Their  properties  are  similar  to  properties  of  ordinary  differential 
equations  with  a lagging  argument  [17]. 

This  set  has  the  only  stationary  point  a=  1 , / = 1 , C = 0.5.  It  corresponds  to  the  steady  supercavity 
(29).  We  linearize  the  equations  in  the  stationary  point  neighbourhood,  representing  the  unknown 
functions  in  the  form: 


a(t)  = l + eae^ , l(t)  = \ + £beM> , C(t)  = — + eceM‘ 

and  saving  the  only  terms  of  order  O(e) . Here,  u = X + jk,  k = col0  / is  the  reduced  frequency  of  the 
oscillations;  a,  b,  c are  constants.  In  this  case  the  integrals  in  the  equations  are  expressed  by  the  modified 
Bessel  functions  of  the  1st  kind  I0(/u/ 2)  and  /,  (/j  / 2)  [27].  As  a result,  we  obtain  the  set  of  the  linear 
uniform  equations  with  respect  to  unknowns  a,  b,  c.  We  do  not  present  it  here  because  its  awkwardness. 
Equating  the  determinant  of  the  equation  set  to  zero,  we  obtain  the  characteristic  equation: 

A [n 2 (/ 2 - 1\ ) + 2ju/0 (21 0 + /, ) + 21 2 ] + \MA  In  2 + 4(j8  - 1)][2/0/1  +M(I0+  /,  )(/0  + 31, )]  = 0,  (44) 

A = 8[mQ0  + y(P-1)]/x. 

where  Q0  = naQll  / 8 is  the  volume  of  the  stationary  pure"  cavity  (29). 

We  suppose  firstly  that  the  mass  of  gas  in  the  cavity  is  constant,  i.e.  y=  0.  In  this  case  the  Eq.  (44) 
contains  the  only  physical  parameter  /) . It  is  easy  to  make  sure  that  )i  =0  is  the  only  real  root  of  the 
Eq.  (44).  Assuming  X = 0 and  using  the  relations: 


where  J0,J , are  the  Bessel  functions  of  the  1st  kind,  we  establish  that  the  Eq.  (44)  has  two  series  of 
alternating  pure  imaginary  roots,  i.e.  frequencies  of  the  fundamental  oscillations  k„  ,k„  ,...  They 
satisfy  the  equations: 


J , 


V 2 7 


= 0, 


/c/1(J02+J12)  + J0(J02-J12)  = 0. 


Corresponding  values  of  the  parameter  are: 

/tf}=  1 + 


(1)_,^2£>0(31n2  \)k2'  Jg(2)=1  + 


22o 


3 n 


n 


Jx(k!  2) 


(45) 


A numerical  analysis  of  the  complex  roots  of  the  Eq.  (44)  shows  that  we  have  X > 0 when 
k?  < k < /c'.2 1 , n = 1,2, ... , i.e.  oscillation  with  such  frequencies  is  unstable.  A graph  of  the  dependence 
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of  X(/3)  is  presented  in  Fig.  26.  A comparison  with  Fig.  11  shows  that  the  distributions  of  the 
characteristic  equation  roots  is  qualitatively  similar  for  the  dynamic  models  of  both  the  two-dimensional 
and  the  axisymmetric  gas-filled  cavities. 


Fig.  26.  Real  part  of  the  characteristic  Eq.  (44)  roots 


When  1 < p < /i,'1’ , the  steady  cavity  is  asymptotically  stable,  and  when  /J  > /31(1) , it  is  unstable. 
When  the  parameter  /I  increases  and  passes  through  the  values  /3*1’ , the  oscillations  with  frequencies 
are  excited.  We  may  suppose  that  if  several  frequencies  simultaneously  exist  for  the  given  value  of 
[5  , the  frequency  corresponding  to  the  highest  coefficient  of  growth  X becomes  dominate  when  the 
nonlinear  self-induced  oscillation  of  the  cavity  is  establishing.  This  allows  to  estimate  the  limits  of  modes 
of  the  cavity  oscillation.  A graph  of  dependence  k(  (i)  is  presented  for  the  first  five  modes  in  Fig.  27 
(solid  lines). 

When  y 0 as  in  the  case  of  axisymmetric  cavities,  number  of  the  fundamental  cavity  frequencies 
becomes  finite.  When  y > 0.2153,  the  cavity  is  asymptotically  stable  for  any  values  of  (5  . According  to 
the  data  [6]  for  the  stationary  cavities  past  a flat  plate  perpendicular  to  the  mainstream  y = 2.5  • 10 -4  and 
effect  of  taking  into  account  the  gas-leakage  variability  is  small. 


Table  1.  Comparison  with  experiment  (the  Is'  M. Tulin's  scheme) 


Unbounded  flow 

Jet  H = 0.254  m 

Experiment  [6] 

W 

W 

W 

lj/H 

/,//„ 

iC 

W 

Pn 

kN 

3.081 

15.40 

4.810 

7.088 

1.875 

0.589 

2.226 

9.485 

5.06 

6.15 

15.40 

42.70 

11.28 

13.38 

2.125 

0.515 

8.418 

22.48 

16.1 

12.9 

42.70 

83.60 

17.74 

19.64 

2.000 

0.551 

23.97 

46.51 

33.5 

18.8 

83.60 

138.1 

24.12 

25.90 

2.000 

0.551 

46.51 

76.55 

55.5 

22.9 

138.1 

206.5 

30.47 

32.17 

1.875 

0.589 

81.74 

122.0 

88.7 

30.3 

The  experiments  [6]  were  carried  out  in  a falling  jet  with  the  width  H = 25.4  cm  at  the  decreased 
external  pressure  . Closeness  of  the  free  boundaries  of  the  flow  results  in  decreasing  both  the  length 
1 1 and  the  volume  Qj  of  the  steady  cavity.  According  to  the  relations  (45)  this  should  affect  on  the 
parameter  (5  values.  Using  the  results  of  the  work  [28],  we  obtain  the  estimation  of  influence  of  the  jet 
boundaries  on  the  "pure"  supercavity  volume  when  a — > 0: 


2l~Il 

Q0  /„  ’ 


H 


sinh  — 


H 
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Table  1 represents  theoretical  values  of  the  parameters  /3  and  k for  unbounded  flow  which 
correspond  to  the  limits  of  the  cavity  pulsation  modes  The  parameter  /3  values  are  recalculated  by 
formulae  (45)  with  taking  account  of  influence  of  the  jet  boundaries. 

A graph  of  the  theoretical  dependence  k(  /3 ) with  taking  account  of  the  jet  boundaries  for  five  modes  of 
the  cavity  pulsation  are  plotted  in  Fig.  27  by  dotted  lines.  Experimental  points  from  [6]  also  are  plotted 
there  by  circles.  We  have  a good  agreement  with  taking  account  of  the  experimental  data  dispersion. 

k 
30 

£5 

£0 

15 

10 

0 100  200  p 

Fig.  27.  Modes  of  the  two-dimensional  supercavity  pulsation 

In  our  work  [21],  we  have  investigated  stability  of  the  2-D  ventilated  supercavities  in  a free  jet  with 
using  the  2"a  linearized  M. Tulin's  scheme  with  the  infinite  wake  and  limited  pressure  at  the  infinity.  We 
shown  that  the  dynamic  properties  of  both  the  cavitation  schemes  are  identical. 

6.3  The  2-D  supercavity  shape  at  sinusoidal  time  perturbations 

Now  we  describe  a method  of  calculation  of  the  length  and  the  shape  of  the  unsteady  2-D  supercavity 
with  using  the  integral  Eqs.  (36)  and  (39). 

The  set  of  the  equations  (36),  (39)  and  (4)  is  nonlinear  one,  if  the  cavity  length  l(t)  is  considered  as  an 
unknown  time  function.  The  cause  is  that  variation  of  the  function  l(t ) has  order  lit )~  0(\ ) even  if 
<7  ~ 0(e)  and  (pv(x,t ) ~ O(e) , where  e is  the  small  parameter. 

We  consider  a practically  important  case  of  periodical  dependence  of  the  flow  on  time  [22].  Let  the 
upper  wedge  side  oscillates  and  deforms  sinusoidally 

f(xJ)  = ax  + k Re  {/*  ( x)e'kt } , (46) 

where  a~  k ~ 0(e) . Here  and  below,  values  which  are  complex  with  respect  to  j,  are  marked  by  the  star. 
These  are  named  the  complex  amplitudes. 

In  the  Eq.  (36),  we  replace  the  time  derivative  inside  the  integral  and  use  exponential  representation 
for  the  time  functions  with  variation  of  order  O(e) : 

q(x,t)  = aq0(x,/)  + tcRe{q  (x)e’kt },  o(t)  = acr0(/)  + vRe{crV/fe }, 

F(x , t)  = ccF0  (x, I)  + k Rc{F  (x)e^kt }. 


(47) 
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The  first  terms  in  the  relations  (47)  represent  quasi-stationary  components  of  the  solution,  which  depend 
on  time  only  via  l{t).  The  second  terms  represent  non-stationary  perturbed  components  of  the  solution. 

One  can  show  that  replacing  the  time  derivative  inside  the  integral  and  also  representation  of  the 
solution  in  the  form  (47)  is  equivalent  to  neglecting  some  terms.  The  effect  of  the  latter  is  essentially 
decreased  with  increasing  the  average  cavity  length  l. 

Substituting  the  expressions  (47)  into  the  Eqs.  (36)  and  (39),  we  obtain  a set  of  the  singular  integral 
equations  with  respect  to  the  complex  amplitudes  q*(x)  = q^(x)  + jq2(x ) and  a*  = <7,  + j<J2 : 


+ jk  In  | x - s 


s-na 


■■  A*  (x) , 


id) 

JV(s)cfe  = A*2, 


(48) 


(49) 


where  A*(x) , A2(x)  are  known.  When  k = 0,  we  obtain  from  the  Eqs.  (48)  and  (49)  the  equations  to 
determine  the  quasi-stationary  components  q0(x),  a0. 

We  consider  the  cavity  length  l(t)  as  a free  time  parameter  and  determine  it  at  the  sequential  times 
t(n)  = t(n~l)  + At  by  numerical  solving  the  equation  of  the  mass  of  gas  in  the  cavity  balance  (4).  In  this 
case  the  unsteady  cavity  volume  Q(t{n))  is  calculated  by  integrating  the  expression  (38).  Using  the  Eq. 
(38)  for  the  sinusoidal  oscillation,  we  have 

X 

F\x)  =^e~jkx  J q\s)ejhds  . (50) 

o 

For  each  iteration,  both  the  quasi-stationary  and  the  non-stationary  parts  of  the  solution  are  calculated  at 
the  same  fixed  value  l(n)  from  the  set  of  Eqs.  (48)  and  (49)  by  the  numerical  method  of  discrete 
singularities  [24]. 


z,  X 

H= ► 

1 z=Vx 

x - sources 
o - control  milestones 


Fig.  28.  Discretization  of  the  problem 


The  numerical  method  of  discrete  singularities  consists  in  approximation  of  the  integral  Eqs.  (48)  and 
(49)  by  a set  of  linear  algebraic  equations.  In  this  case  the  continuous  distribution  of  sources  along  the  x- 
axis  is  replaced  by  discrete  one,  and  then  the  quadrature  formula  of  rectangles  is  applied.  Change  of 
variables  „t— >z2,  s— >A.2  is  preliminarily  realized  in  the  integrals  to  improve  the  method  convergence 
[24].  The  cavity  projection  is  divided  into  M equal  intervals.  A point  source  and  a collocation  point 
(control  milestones)  are  located  in  each  interval.  The  boundary  condition  (24)  is  satisfied  at  the 
collocation  points.  Order  of  location  of  singularities  and  collocation  points  z,  is  determined  by  the 
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class  of  the  solution  of  the  singular  integral  Eq.  (48).  This  solution  must  be  limited  at  the  point  of  the 
cavity  separation  x = 1 , and  it  is  unlimited  at  the  point  x = /(n) : 

JJw  _ i 

z,  =l  + 4z(/-0.75),  Xj  =1  + Az(j  -0.25)  Az=— , i,  j = l,  2,...,M.  (51) 

J M 

A scheme  of  the  problem  discretization  is  shown  in  Fig.  28  when  M=  3. 

As  a result,  separating  the  real  part  and  the  imaginary  part,  we  obtain  a set  of  2 (M  + 1)  linear  algebraic 
equations: 


M 

7-1 


M 

7=1 


in) 

- A- 

( „<"> 


' 0 ' -i  >' 


I r,  a2  “ 7 ,2 


( = 1.2 W, 


( = 1,2 (W, 


Fig.  29  shows  the  characteristic  wave-like  shape  of  the  cavity  past  the  wedge  when  its  sides  are 
performing  pitching  oscillations  about  the  nose 

f(- v)  = -r 


for  three  values  of  the  reduced  frequency  k. 


Fig.  29.  The  cavity  shape  in  dependence  on  the  oscillation  frequency:  /0  = 10.0,  K =0.1 

In  each  case  the  cavity  shape  was  calculated  by  the  Eqs.  (47)  and  (50)  at  the  time  tk  when  Htk)  = l0 
for  convenience  of  comparison.  A shape  of  the  closed  stationary  supercavity  when  k = 0 is  shown  by 
dashed  line.  When  k>  0,  the  cavity  is  unclosed. 

Character  of  the  cavity  deformation  is  the  same  for  different  types  of  oscillations  of  the  wedge. 
Kinematic  waves,  which  are  formed  by  oscillation  of  the  points  of  the  cavity  boundary  separation,  move 
along  the  cavity  with  the  velocity  V „ . Their  amplitude  increases  approximately  by  linear  law. 
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Graphs  of  the  functions  I(t ) and  Q{t),  which  are  calculated  for  only  period  of  the  oscillation  and  various 
values  of  the  relative  amplitude  k = k/cc,  are  given  in  Fig.  30.  The  wedge  sides  are  wave-likely 
deformed  according  to  the  law  (for  the  upper  side): 

f{x)  = xe~jkx . 

One  can  see  that  the  cavity  length  oscillation  distinguishes  more  and  more  from  sinusoidal  one  when  the 
amplitude  K increases.  The  same  occurs  when  the  reduced  frequency  k increases.  The  functions  1(1)  and 
Q(t)  become  discontinuous  when  k and/or  k exceed  the  some  critical  values. 

Fig.  31  shows  time  dependence  of  the  axisymmetric  cavity  length  and  volume  that  calculated  with  the 
code  PULSE.  Flere,  a cause  of  the  cavity  unsteady  perturbation  is  the  ambient  pressure  oscillation: 

P(t)  = P0(1  + k coskt) . 


Fig.  30.  Influence  of  the  oscillation  amplitude  on  the  functions  l(t)  (a)  and  Q(t)  (b):  l0  = 6.0,  k = 0.5 

(wave-like  deformation  of  the  wedge) 


Fig.  31.  Flistory  of  the  length  (a)  and  the  volume  (b)  of  the  axisymmetrical  supercavity:  f30  = 1.0;  k=  2.5 

(the  ambient  pressure  pulsation) 

A comparison  of  graphs  in  Figs.  30  and  31  shows  that  shape  of  the  graphs  of  l(t)  and  Q(t)  for  the 
axisymmetric  and  the  two-dimensional  cavities  is  qualitatively  similar  in  the  case  of  nonlinear 
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oscillations.  When  the  oscillation  frequency  and  amplitude  increases,  their  shape  more  and  more  deviates 
from  sinusoidal  one  having  a shape  of  a "falling  wave".  When  the  frequency  and/or  the  amplitude  exceed 
some  critical  values  kcr  and  Kcr , the  functions  /(/)  and  Q(t)  become  discontinuous. 

The  same  kind  of  behaviour  of  the  unsteady  supercavities  is  observed  at  the  experiments  for  both  the 
axisymmetric  and  two-dimensional  cavities  [29], 


References 

1.  Savchenko,  Yu.N.,  Vlasenko,  Yu.D.  and  Semenenko,  V.N.,  Experimental  Study  of  High-Speed 
Cavitated  Flows,  Int.  J.  of  Fluid  Mechanics  Research , 1999,  26,  No.  3,  pp.  365-374. 

2.  Savchenko,  Yu.N.,  Semenenko,  V.N.,  and  Putilin,  S.I.,  Unsteady  Supercavitated  Motion  of  Bodies, 
Int.  J.  of  Fluid  Mechanics  Research,  2000,  27,  No.  1,  pp.  109-1 37. 

3.  Logvinovich  G.V.,  Hydromechanics  of  Free-Boundary  Flows,  Halsted  Press,  1973. 

4.  Epshtein,  L.A.,  Methods  of  Theory  of  Dimensionality  and  Similarity  in  Problems  of  Ship 
Hydromechanics,  Sudostroenie  Publishing.  House,  Leningrad,  1970  [in  Russian]. 

5.  Korolyev,  V.I.,  Vlasenko,  Yu.D.,  and  Boyko,  V.T.,  Experimental  Investigation  of  the  Cavity 
Development  at  Unsteady  Gas  Cavitation,  Gidromehanika,  1973,  No.  24,  pp.  79-83  [in  Russian], 

6.  Silberman,  E.,  Song,  C.S.  Instability  of  Ventilated  Cavities  , J.  of  Ship  Res.,  1961,  5,  No.  1,  pp.  13- 
33. 

7.  Michel,  J.M.,  Ventilated  Cavities.  A Contribution  to  the  Study  of  Pulsation  Mechanism,  Unsteady 
water  flows  with  high  velocities,  Proc.  of  International  Symp.  IUTAM,  Nauka  Publishing  House,  1973, 
pp.  343-360. 

8.  Paryshev,  E.V.,  Theoretical  Investigation  of  Stability  and  Pulsations  of  Axially  Symmetric  Cavities, 
Trudy  TsAGI,  1978,  No.  1907,  pp.  17-40  [in  Russian]. 

9.  Paryshev,  E.V.,  Numerical  Modeling  of  Ventilated  Cavity  Pulsations,  Trudy  TsAGI,  1985,  No.  2272, 
pp.  19-28  [in  Russian]. 

10.  Knapp,  R.T.,  Daily,  J.W.,  Hammitt,  F.G.,  Cavitation,  McGraw-Hill,  New  York,  1970. 

11.  Zhuravlyev,  Yu.F.,  Methods  of  Theory  of  Perurbations  in  Three-Dimensional  Jet  Flows,  Trudy 
TsAGI,  1973,  No.  1532,  pp.  8-12  [in  Russian]. 

12.  Savchenko,  Yu.N.,  Semenenko,  B.H.  Wave  Formation  on  the  Boundaries  of  Cavities  Forming  at 
Water  Entry  of  a Disk  and  Cones,  In:  Problems  of  High-speed  Fluid  Mechanics,  Chyvashi  University 
Press,  Cheboksary,  1993,  pp.  231-239  [in  Russian]. 

13.  Logvinovich,  G.V.,  Syeryebryakov,  V.V.,  On  Methods  of  Calculations  of  Slender  Asisymmetric 
Cavities,  Gidromehanika,  1975,  No.  32,  pp.  47-54  [in  Russian]. 

14.  Semenenko,  V.N.,  Computer  Simulation  of  Dynamics  of  Supercavitating  Bodies,  Prykladna 
gidromehanika,  2000,  2,  No.  1,  pp.  64-69  [in  Russian]. 

15.  Semenenko,  V.N.,  Computer  Modeling  of  Pulsations  of  Ventilated  Supercavities,  Int.  J.  of  Fluid 
Mechanics  Research,  1996,  23,  Nos.  3 & 4,  pp.  302-312. 

16.  Semenenko,  V.N.,  Computer  Simulation  of  the  Unsteady  Supercavitating  Flows,  High  Speed  Body 
Motion  in  Water  (AGARD  Report  827),  Proc.  Fluid  Dynamics  Panel  Workshop,  Kiev,  1997. 

17.  Elsgolts,  L.E.,  Norkin,  S.B.,  Introduction  in  Theory  of  Differential-Delay  Equations,  Nauka 
Publishing  House,  Moscow,  1971  [in  Russian]. 


12-30 


18.  Logvinovich,  G.V.,  Yakimov  Yu.L.,  High-Speed  Immersion  of  Bodies  in  Fluid,  Unsteady  water 
flows  with  high  velocities,  Proc.  of  International  Symp.  IUTAM,  Nauka  Publishing  House,  Moscow, 

' 1973,  pp.  85-92. 

19.  Semenenko,  V.N.,  Instability  of  a Plane  Ventilated  Supercavity  in  an  Infinite  Stream,  Int.  J.  of  Fluid 
Mechanics  Research,  1996,23,  Nos.  1 & 2,  pp.  134-143. 

20.  Semenenko,  V.N.,  Instability  and  Oscillation  of  Gas-Filled  Supercavities,  Proc.  Third  International 
Symp.  on  Cavitation,  Grenoble  (France),  1998,  2,  pp.  25-30. 

21.  Semenenko,  V.N.,  Instability  of  a Plane  Ventilated  Supercavity  in  an  Free  Jet,  Prykladna 
gidromehanika,  1999,  1,  No.  2,  pp.  45-52  [in  Russian], 

22.  Semenenko,  V.N.,  Calculation  of  Shape  of  Plane  Supercavities  at  Harmonic  Oscillations,  Prykladna 
gidromehanika,  2000,  2,  No.  3,  pp.  87-93  [in  Russian], 

23.  Tulin,  M.P.,  Supercavitating  Flows  - Small  Perturbation  Theory,  J.  of  Ship  Research,  1964,  7,  No.  3, 
pp.  16-37. 

24.  Yefremov,  LI.  Linearized  Theory  of  Cavitation  Flow,  Naukova  Dumka  Publishing  House,  Kiev, 
1974  [in  Ukrainian], 

25.  Gakhov,  F.D.,  Boundary  value  problems,  3rd  Revised  and  Supplemented  Edition,  Fizmatgiz 
Publishing  House,  1977  [in  Russian], 

26.  Semenenko,  V.N.,  Unsteady  Flow  Calculations  Past  Ventilated  Hydrofoils,  Proc.  Seventh 
International  Conf.  on  Numerical  Ship  Hydrodynamics,  Nantes  (France),  1999. 

27.  Abramovitz,  M.,  Stegun,  I.  (Editors),  Handbook  of  mathematical  functions  with  formulae,  graphs  and 
mathematical  tables.  National  Bureau  of  Standards,  Applied  Mathematics  Series,  1964. 

28.  Epshtein,  L.A.,  Lapin,  V.M.  Approximate  calculation  of  influence  of  Flow  Boundaries  in  Two- 
Dimensional  Problem  and  past  an  Axisymmetric  Body,  Trudy  TsAGI,  1980,  No.  2060,  pp.  3-24  [in 
Russian]. 

29.  Nishiyama,  T.,  Unsteady  Cavity  Flow  Model  for  Two  Dimensional  Super-cavitating  Hydrofoils  in 
Oscillation,  Technology  Reports,  Tohoku  Univ.,  1982,  46,  No.  2,  pp.  199-216. 


